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Abstract 

Recent analyses suggest that in TeV scales that will be made accessible at the LHC copious 
amounts of color scalar parton bound states may be produced. Would this be the case, 
the scalars would leave long enough to interact and this could lead to new physics. These 
interaction could be direct, but also take place with a mediation of the dense parton medium 
through which they propagate. Since multiple processes would have to be included, the latter 
case is too convoluted to be treated with perturbative methods applied to the Standard Model. 
We explore a possibility of describing such interaction by a nonlocal Lagrangian which is an 
infinite polynomial in the field variables, momentum and mass. We treat all scalars as identical, 
use a 0(Af)-symmetric Lagrangian, where N is the number of scalars and discuss the problem 
in expansion. Nonrenormalizable by all traditional criteria, such model still requires 

only a limited set of multiplicative renormalizations, provided that the parameters in the 
Lagrangian are not treated as an infinite set of independent coupling constants, but as finite 
expansion coefficients of the Lagrangian in powers of field variables, mass and momenta. The 
necessary constraints on the relative magnitudes of these coefficients can be determined order 
by order in a double series in the single coupling constant and 1/7V. 
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1 Introduction 



The idea of extending the Minimal Standard Model to include higher symmetries and richer 
particle spectrum is almost as old as MSM itself. For example the proposition of a four color 
symmetry in which lepton numbers are considered as the fourth color dates to 1974 [IJ. The 
major obstacle in increasing the numbers of particles in the model was however related to the 
fact that this generally leads to unacceptably strong flavor-changing neutral currents. This flaw 
turned out to be absent in the case of minimal flavor violation [2]. It effectively reduces the 
permitted extension of MSM to fields that couple to partons as scalars that carry gauge or color 
quantum numbers. The number of representations permitted by the requirement that scalars 
Yukawa-couple to fundamental fermions is very limited, and for SU{3)c x SU{2)l they include 
(8, 2), (6, 3), (6,1), (3, 3), (3, 2), (3, 2), (3,1), (1,3), (1,2), and (1,1) representations. Of these, even 
though present in the liteature from the era of technicolor [3], the color octet is a particularly 
intriguing candidate, as it has the same weak quantum numbers as the Higgs doublet [4], [5]. At 
LHC energies, the total cross section for production of such such scalars, with masses ^ l.iTeV is 
large enough for them to be observed in large quantities as heavy resonances resulting from parton 
merging. [6] 

The color scalar bound states are expected to live log enough to ineract. It was even suggested 
that they can produce their own (color singlet) bound states, [7j, thus there is no fundamental 
argument against the conjecture that scalars from all mentioned representations significantly co- 
participate in partonic interactions. For example, it has been shown, on phenomenological grounds, 
that the decays of the color sextet scalars could be responsible for a significant contribution to the 
production of top quark pairs even below the 1 TeV region of scalar mass [S]. 

Among the allowed vertices involving composite scalars, besides the obvious qqs and ggs cou- 
plings we also find the ggss vertex. This implies the possibility of building purely scalar loops, 
and also the production of cascades of gluon pairs by energetic color scalars propagating in the 
parton medium. The latter could supplement (but not compete with) the known mechanisms of 
production of gluon cascades by sequential emissions of single gluons by energetic pomerons. 

A separate question is that of the nature of mutual interaction between colored scalars. The 
simplest proposition would be to use the contact interaction potential term in the form of 

with the summation carried over all scalars (from all relevant, rather than restricted to one chosen 
representation) . The weakness of such proposition is that it does not involve indirect scalar-scalar 
multiple interactions that take place with the mediation of the dense parton medium. Because 
of the sheer number of possible processes the problem is too complex to be treated perturba- 
tively, however a nonlocal effective interaction might provide an approximate description of the 
phenomena involved. 

The goal of this contribution is to seek a possible form of such effective interaction. We will 
employ a.l/N expansion (with N representing the total number of composite scalars, as opposed 
to the dimensionality of the group representation). Further, as the first approximation we shall 
ignore individual differences between the scalars (charge, mass, etc), which are secondary for the 
problem, and consider a nonlocal, 0(7V)-symmetric interaction of identical neutral fields. 

Our goal is to explore infinite-polynomial interactions of scalar fields (that also include higher- 
order derivative terms), without treating each coefficient of the expansion as an independent in- 
teraction coupling that requires a separate renormalization constant. Instead, we will treat them 
as finite Taylor expansion coefficients (in the field, mass and momentum variables) of the total 
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interaction Lagrangian. In no case these coefficients sliould be regarded as expansion parameters 
of a predetermined form-factor that assures finiteness of the amphtudes. The procedure will be 
carried order-by-order in 1/A^ and in compliance with the standard power counting rules in the 
coupling constant. 

The form of the suggested Lagrangian is typical for nonrenormalizable interactions. However, 
since we will be working under a crude assumption that all scalars from several possible represen- 
tations interact identically, we can not assume that this effective interaction could possibly be valid 
beyond several orders in where the approximate 0{N) symmetry remains valid. Therefore 

we will seek an exact form from of the interaction Lagrangian (in terms of its expansion coefficients 
in all relevant field, mass and momentum variables) from the requirement that, order by order in 
1/N, and m^, it is rcnormalizable in the sense that it requires only a finite number of standard 
renormalization constants, similar to that of the $^ model. 

We will use the dimensional renormalization technique and will demonstrate that it is possible 
to algebraically determine the necessary relations between the expansion coefficients that lead to 
mutual cancellations between the harmless (polynomial) poles. We also offer an argument that our 
procedure does not lead to harmful logarithmic poles. The result is almost unique, in the leading 
order only one of the coefficients of expansion remains free. In this sense, the condition that the 
number of renormalization constant remains finite defines the form of the integration Lagrangian. 
Still, we do not make a claim that the method presented here is much more than a tool for effective 
models. 

An alert reader may be alarmed that in our procedures the order of separating pole parts and 
taking the large- iV limit seems to be reversed. Indeed, in the case of cutoff procedures applied to 
U{N) models, the limit of large N is taken at a finite value of the cutoff parameter. If the infinite 
limit of the cutoff parameter is taken first, the procedure is generally invalidated at any finite 
(although large) order in N. In the case under consideration the order of operations is immaterial. 

The singularities produced by the dimensional method applied to the proposed scalar 0{N) 
interaction are simple poles with the residues coming from a limited number of diagrams which is 
much smaller than N. All but one of the summations involved are finite. The only exception is the 
well discussed and accommodated 'bubble chain' of graphs, but this contribution is of the order 
of 0{N'^). In higher orders in the procedure is therefore perturbative. Indeed, the numbers 
of contributing graphs increase because of the contributions from diagrams topologically identical 
to the lower-order ones but with redirected A'^- number flow. However, this effect is algebraic and 
is not affected by the change of the order of summation and renormalization. Further, as just 
explained, the proposed model itself is not expected to retain its validity in such orders of the 
expansion. 

The question of reconciling the order of resummation and renormalization in the 0{N)^^ model 
within the cutoff, rather than dimensional, method has already been more formally discussed 
leading to similar conclusions. [9] 

The choice of the 1/N expansion is not accidental; our task would not be possible within 
confines of the conventional perturbation theory which in every order of the coupling contant 
mixes vertices with arbitrary numbers of external legs and any possible momentum dependence. 
In contradistinction, the 1/N power counting rules establish unique relations between the order of 
approximation and the maximum number of external legs in diagrams. The infinite summations 
of graphs in the lowest orders of the expansion provide a nonperturbative ground state which 
serves as an alternative for the conventional free-field solution. Last but not least, the powers of 
1/A^ provide us with a finite, unaffected by renormalization procedures, non-running expansion 
parameter for order-by-order adjustment of the required magnitudes of the expansion coefficients 
in the Lagrangian. 
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We end this introduction with a few remarks on nonrenormahzable interactions which are 
beheved to be able to provide Hmited but useful insight into the nonperturbative realm outside 
the limits of the asymptotic freedom. 

Already a few years after the discovery of renormalizability of non-abelian theories, Weinberg 
argued [lOj that, despite the known problems with infinite number of renormalization conditions, 
with judicious approach, some nonrenormahzable models may still hold predictive powers if the 
nonrenormalizability is understood as 'a necessary evil', in the sense that the solutions to a number 
of initial terms of the expansion of the effective model mimic the behavior of underlying renor- 
malizable and consistent theories. More recently, the argument was made that it might be 
possible 'not to restore renormalizability in the Dyson (dimensional power counting) sense, but to 
live without it' . With the imposition of proper a priori constraints on the bare action, a gauge 
theory that is nonrenormahzable in the usual sense, can nevertheless be renormalized by introduc- 
tion of an infinite number of renormalization terms in the quantum action while maintaining its 
predictive powers. 

In a different context of nonrenormalizability of scalar interactions in more than four dimen- 
sions, Klauder [1^ (see also his previous work in Ref. [I3|), argued that nonrenormahzable inter- 
actions are in fact discontinuous perturbations, and therefore adding counterterms to the original 
free theory may be an irrelevant procedure. Instead, he proposed evaluating solutions about an 
alternative pseudo-free theory. 



2 The model 



The Lagrangian density of our model is a generic 0(A'')-symmetric infinite-polynomial interaction 
of A'' real (easily generalizable to complex) scalar fields 

Ca^^^d^'^'^d,^'' -^^ml^''^'' - Cs- Ci, (1) 

where Cs (for future convenience written as an interaction term) represents the nonlocal part of 
the propagator, 

k+l=2 

In addition to /i, which is the conventional mass unit scale, Lq in ^ represents the scale of 
nonlocality range. The interaction Lagrangian has the form 

where 



D^^^d'i^Cd^^.timlfiLlt^K (4) 

Since ^ exhibits all marks of a nonrenormahzable theory, some clarifications are necessary 
before we continue. The renormalizable ■^(<f>''<I''')| model has three types of divergent operators 
and three renormalization constants, Zg, Z^, and Z$. In the established notation. 
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1+ E ^"-^^j.^;^ 



(6) 



and 



1 

In the model under consideration, the parameters a,yj, 6,^^, and Cjyj are also implicit series in 
1/iV. For sake of transparency of the notation we will not introduce indices of the expansion in 
powers of this parameter. 

The role of derivative terms in the Lagrangian is not to a priori introduce functions of mo- 
menta that modify the vertex or propagator functions so that the standard power counting rules 
no longer apply and the divergence problem softens or disappears, as in the early 'finite theory 
without renormalization' concept of nonlocal interactions. Instead we will start from the working 
assumption that in the initial order in the asymptotic behavior of the propagators is canonical, 
and that the local-like power counting rules for the vertices apply in the asymptotic domain. Then 
we will show by construction that it is possible to to introduce constraints on the expansion pa- 
rameters so that thanks to mutual cancellations there is no need for the proliferation of the number 
of renormalization constants. In a sense, the method resembles quantum mechanical procedures 
in which the physical solutions are found by choosing the expansion parameters of a wave function 
on the requirement of its normalizability. We do not make any claims about consistency of the 
resummed "full theory". While the reality may turn out to be more generous, we seek only a tool 
for an effective model that behaves as a renormalizable theory through a finite number of orders 
of the 1/N expansion. 

The computational techniques that we will use are not new; the 1 /iV-expansion originates from 
the seventies, |14| . |15| . the infinite-polynomial Lagrangians (without the derivative terms) were 
first extensively explored by Schnitzer, [T^ and later also by this author, who applied the method to 
the non-leading order and bounded (in the field variable) interactions [TH] , [T^ . The regularization 
of choice will be dimensional, [20] even though we will not attempt to provide a proof that all 
conditions of consistency are met to all orders. However we will offer an argument that, due to 
infinite summations, the procedure is not affected by the problem of logarithmic poles due to the 
the graph topology which includes the 0((1/7V)°) 'bubble-chain' contributions in every diagram 
category. 



3 The leading order 

The principal \/N power-counting rule for polynomial 0{N) scalar Lagrangians is that every 
vertex with n pairs of external legs contributes a factor of (1/A^)"~^, and every closed loop that 
involves the summation over the field indices contributes a factor of A^.[Tn] One consequence of 
these rules is that no graph with 2n external legs can appear in the {\/N)^~^ order. 
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This property of the \/N expansion is key for our purpose since, unhke in the conventional 
perturbation expansion, each subsequent order adds only one new class of divergent amplitudes 
(with two more external legs than the maximum allowed in the preceding order). Accordingly, the 
dominating order contains only the vacuum graphs which we shall ignore here, even though they 
could be important in cosmological applications. The leading order (LO) contains only the self- 
energy functions, the next-to leading order (NLO) additional propagator graphs plus four-point 
functions, and so forth. 




Figure 1: "Petal" graphs that contribute to the effective four-point vertex. All vertices in the 
following graphs are effective vertices, and the heavy dot is not used. 



Still, vertices of order higher than ($^)^ do contribute to the LO, but only as parts of the petal 
graph series presented in Fig. 1 which form the 'effective' ($^)^ vertex. The petal loops do not 
carry external momenta, and they do not produce overlapping divergences. Therefore, despite the 
presence of derivative interaction terms, we will be able to sum them directly. In purely polynomial 
models without derivative couplings, such summation is equivalent to taking a derivative of the 
interaction potential with respect to at a certain value of the effective field. In the case under 
consideration we must proceed in a more deliberate manner. 




Figure 2: Daisy graph contributions to the propagator. 



Consider the simplest single- loop propagator diagram shown as the first graph in Fig 2. It 
actually represents an infinite sum of graphs that include vertices with all possible powers of mass 
and momenta with their respective coefficients d\i . The success of the proposed procedure depends 
on whether it is possible to remove all the divergences of these graphs by a single renormalization 
of mass as in (jni). 

Since we are interested only in the general renormalization properties and work on the assump- 
tion that the ordinary power counting rules apply in the asymptotics, we will ignore the nonlocal 
part of the 2-point function and replace it by an ordinary local canonical propagator. Then the 
first graph of Fig. 2 is a sum of expressions of the form 

With help of the identity r(fc -I- f )r(l - A: - f ) = (-l)'=r(f )r(l - f ), and the fact that so far 
5o — 9r, we find that Eq.® yields 

4.-2^(^l-2j(^-o(^^j fcHiK^o) ■ (9) 



6 



Since the last expression involves arbitrary powers of mass with their respective coefficient 
factors, there are alternative ways of dealing with the divergence contained in ([9]). One possibility 
is to use the fact that r(l - f ) = + 7 - 1 + 0{n - 4), take 



,2 _ / 1 "00 9R 



(10) 

and simultaneously introduce a renormalization constant for the interaction range Lq so that 
ttiqLq = m\L\. 

The other option is to retain the condition (jlOp . but instead of renormalizing Lq, compensate 
the effects of higher powers of ttiq by imposing constraints on the finite coefficients d^,""* . While the 
latter approach is more compatible with the proposed technique, it leads to conditions that link 
d^^i to yet unexplored d!]^i . We relegate the explanation of the details of this fact to the Appendix. 

In any case, the necessary constraints on d^"'' for n > 1 can be determined solely from the 
assumption that the asymptotic behavior of the nonlocal propagator is canonical, without pos- 
sessing any detailed knowledge of its exact form outside the asymptotic region. The procedure is 
conditional on whether or not the recurrence conditions for the coefficients that define the series 
expansion of the propagator reproduce the canonical behavior at infinity. Would it not be the case, 
then the first procedure would have no alternative. Nonethless, the canonical behavior of the LO 
propagator might still be imposed by hand with an appropriate definition of Cs- 

Hence, we choose to renormalize i, and in analogy to ([6]) we define 

(11) 



1 

ly—lj—i/ 



where in the leading order 



hi = -^11- 

The relation Z^Zm = 1 does not have to hold in higher orders in 1/iV. 

In the following we shall follow as closely as possible the standard renormalization procedures 
developed for renormalizable 1/N models. Accordingly, we shall use the expansion for the 
summation of the dominating classes of graphs while, unless infinite summations of the 0{N'^) 
order are involved, retain the standard perturbation expansion rules for discarding terms of higher 
order in gji. In the case under consideration, the pole parts of all but one Z„i are suppressed by 
Zl, and the Lq associated with momenta are replaced by Lfj since their divergent parts are of 
higher powers in g^. Therefore, for the purpose of LO calculations the Lagrangian effectively takes 
a somewhat simpler form 



^S-^ml Y: ^(^4'l+i$'^(a''5.)'=(m|,)'(L?,,)'=+'$" 



-(a'-cp-a^^^)^ _4l)<i>-(a''a.)'=-l(L2^)'=-lci>^ (12) 



k=2 
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Also 



D^^>^di]\d-^d,r{miyiLj,f+K (13) 

Unlike in a conventionally renormalizablc <i>^ interaction, the mass renormalization alone does 
not suffice to guarantee the removal of infinities. The process involves summations of an infinite 
number of terms that contribute to the residue of the pole, as well as of the finite parts of con- 
tributing diagrams. Further, there are two underlying series summations, one in powers of mass 
and momenta in the vertex at the base of the loop, and another for a series of daisy diagrams 
in Figs. 1 and 2, that make up the effective propagator and (<i)°$'')2 vertex. Both series must 
be summable. Otherwise, one would face infinities of purely algebraic nature which would not be 
curable by standard field renormalization techniques. 

The total pole part of the sum of all contributions to the first graph in Fig. 2 equals 



9r 



7(2) 
^kl 



„_4'"o (4^)2 



2 7-2 \k+l 



k,l 



(14) 



With the simplifying notation 



F«(p^m^)^^|^(p2)'=(m|,)'(L|)'=+^ (15) 
fc,; 

the condition of finiteness of the sum (|14p translates to the condition that T^'^^ (toIj, mj^) is finite. 
With the sum of all vertex terms included, instead of (fTU|) . the coefficient of the first term in Z,„ 
would have to be replaced by 



b'n = j^V('Hml,ml). (16) 

This renormalization guarantees the finiteness of individual petal loops in Fig. 2, but the sum 
of the finite parts of all the daisy graphs 



Si? = 5fl X! 



(J-I)! 



^2 2(47r)2 ^ kl l\ 

h.l 



2 \k+l 



i-1 



(17) 



must also be finite. This is not a restrictive condition, since it suffices that the partial derivative 
of the interaction potential with respect to $2 at — mj^ is finite. 

The LO also includes contributions from the stacked- loop diagrams of Fig. 3. The divergence 
of the top loop is removed by the just completed mass renormalization, while the remaining loops 
require renormalization of the four-point function. Technically, this subject belongs to the NLO, 
but since NLO graphs are of the order (1/A^)°, they can appear as subgraphs in any order. The 
contribution to the LO self-energy is simple since the loops do not carry external momenta. 

Each loop in the series of stacked diagrams contributes 



g = 



ijj. 



9o 



2(27r)" 



^ kl 11 si tl 

k,l,s,t 



2 7-2 



,2 r2 \k+s 



(p2 - ml) 



2\2 ■ 



(18) 
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Then, in consequence of the identity r[^2-f-/c)r(fc+f) = ( f ) T ( 2 - f 

we find that 



^-'2{^\'^) ^\:'2)/^^^k\i\s\t\ f-f ^""^^^^ ■ 

The combined residue, as weU as the finite part of (I19|) converge as a power series provided 
that the sum 



E ^If^I + ^ + ^){mlLlf^''~'^' (20) 

converges. Again, this not a very restrictive condition on the specific shape of the interaction 
potential. It suffices that 



represents an expansion of a regular function that is finite at a; = y = mj^Lj^. Then the combined 
pole parts of (fT9|) can be cancelled solely by the coupling constant renormalization. It is worth 
noticing that the above discussion essentially establishes the sufficient condition for the existence 
of the effective vertex shown in Fig. 1. 
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4 Four- point function in NLO 



The essential technical element of the 1/N expansion is the summation of the infinite series of 
bubble-chain graphs shown in Fig. 5, all of which arc of the order of 0((-^)°). Let us temporarily 
postpone the formal question of the chain summation and first consider the renormalization of 
individual loops that constitute the chain. 




p+k 

Figure 4: A single-loop four-point diagram. 

With derivative couplings, there are additional internal momenta factors in loop integrations, 
but the vertices also provide multiplicative factors of external momenta that flow in- or out- of 
the vertex points. The vertices also contain multiplicative mass factors which combine with those 
produced by integrations over the internal momenta. In the graphs under consideration, each 
internal pair of vertex legs contributes to the numerator of the loop integrand and leads to the 
integrals of the form 

2(27r)" ^ fc! /! s! i! ^'"^^^^-^ J ^ {g^ - mg) [{q + - mlY ^ > 

The integral in ((22)) can be written as 

* / ^>71 21 ^^fdxj^ f 'V'^^p'- . ..p^ni\ (23) 

J {q^ - m'oMq + vY- m^] Jo 



where 



(24) 



In the last formula we have used a simplifying notation M"^ — ttiq — p^x, and P — px. The 
singular, as well as the finite parts of the higher-order (in powers of momenta) of the above integrals 
can be evaluated from the lower-order ones with help of the following recurrence rules for /^', 

7-0 + 1); _ ^ ^^"-1 (or.) 

2(a-l) ^^^^ 

and 

/^e+i) = - 2P^^.,,/(^+i)' + M'Pl (26) 
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In the zero order in m?L'^ and p^L"^ , the contribution is proportional to 



If = -TT-z^r (2 - f ) rf-K - p'^ + (27) 

with a pole part of ^34- In the absence of any other divergences, this pole is removed by the 
standard coupling constant renormalization with 



which we shall assume to work also for the remaining divergent terms, subject to proper restrictions 

(2) 

on the values of the coefficients dj^/ . 

In the next (first) order in the combined powers of m^L^ and p^L^, the four-point function 
collects three types of contributions. One comes from graphs with a factor TOq in one of the two 
vertices, another with a factor of (not subjected to the integration) momenta carried by external 
legs, and lastly, one with an (integrated) single factor of the product of internal momenta within 
the loop. Their combined pole parts expressed in terms of the invariants TOq and {pi ■ P2) equal 



^sl^Ll [ ml (24? +<)) +4(pi ■P2)d^^A (29) 



n - 4 (47r)2 

where we have replaced Lq by Lji and go by which is exact up to 0{g^), 

Despite appearances, there is no need to introduce additional counterterms for the two new 
divergent operators in (P^ . This is because the Lagrangian of interaction ^ already contains 
terms proportional to the invariants m|j and {pi ■ P2) multiplied by go = Zggn which also includes 
the -j^^ factor. If their residues, and also the residues of their respective partners from ([^^ cancel, 
no further subtractions are necessary. The required conditions are 

4oH4oV = 4^0^ 4o\ and 4^1^40^)' = {2d^o^+M^^MoJ, (30) 
which are satisfied under only a mildly restrictive condition 



dg^=i, and 4i' =24f (31) 

The exact value for d^^Q in (|3ip should come as no surprise, since the Lagrangian ([3]) contains 
one free parameter too many {g in addition to the complete set of the djf;'' coefficients). 

In the second order, the term of the type can be produced by graphs that include the 
rfgi^-type terms in both vertices, or one d^^J'^JP^ coupled with a d'^Q partner. The p'^ type requires 
a pair of d^g^'-terms in both vertices, or a djo'' factor in one vertex paired with a d^Q counterpart 
in the other. The latter types contribute also to the m^p^ factor and combine with contributions 
from graphs containing pairs of d^g^ and dgi^ ^r d^^^ and dpo'- before, a proper distinction needs 
to be made between the graphs in which powers of momenta are applied to the internal legs and 
affect the loop integration, or apply to external lines and provide only an overall multiplicative 
factor of external momenta. Then, by matching the pole residues of the m^, m^p^ and p'^ of the 
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sum of all graphs with similar terms from ([3]), one can determine the values ofdj/} , where j + k = 2, 
expressed as functions of similar coefficients with j + k < 2. 

This procedure can be extended to higher orders in mass and momentum, and the fact that the 
value of one seed parameter remains free, guarantees that there exists a degree of diversity in the 
class of possible solutions. However, there are still two possible reasons for the iteration to fail. One 
is that, at some order, the system of equations will contradict itself, and the other that for some 
values of n, k, I, the series of the coefficients d^"'' truncates leading to a simple power of momentum 
behavior. Both possibilities can be eliminated by induction. The argument is straightforward and 
we will present it in a simplified mathematical notation. 

i. Assuming that the n-th order in combined powers of rri^ and is solved, and the values of 
'^^kiji-k) ^'^^ known, in the n + 1 order, the equations for d!^^^j,^^^ have the general form 

j(2) j(2) - ^ ^(2) ,(2) .oo'i 

"OO "fe (ll+l-fc) ~ SlfclJ2fe2"jiA;i"j2fe2- \^^) 



and can be arranged in a hierarchy in decreasing order of the power of momentum k^ . 

(2) 

ii. In the top element of the hierarchy, the coefficient dj^jlj^^-^ g ^PP^ars in the contributing dia- 

grams only once, in the first power, coupled with dpg in the opposite vertex. All the remaining 
contributions of the same combined dimensionality in m? and include coefficients of the lower 
order in s + fc = n, determined in the earlier steps of the recurrence. Given the solutions for the 

(2) 

lower order, the coefficient d(„^]^-)Q is therefore uniquely defined. 

(2) 

iii. Next equation in the hierarchy determines the magnitude of the parameter c?„/. As in the 
previous case, there is only one unknown; the right side of (j32p includes the term proportional to 

(2) (2) 

'^(n+i)o^oo ' ^'^^ remaining parameters have been already determined in the lower orders of 

the iteration. The steps can be repeated down the hierarchy ladder until the process is completed 
at fP'^ 

"0(n+l)- 

(2) (2) . . (2) . . 

iv. Since d^i = 2dj^Q , the procedure still contains one free parameter cJ-^q which sequentially 

(2) 

affects the magnitude of every iteratively generated higher-order parameter c?^; . This degree of 

(2) 

freedom guarantees that it is possible to adjust the magnitude d-^g to prevent the truncation of 
the interaction Lagrangian to a finite polynomial. This is because the number of conditions is 

(2) 

countable, while the unrestricted coefficient d\f^ is real. 

f2) f2) 

It bears emphasis that, given dgj^ , all subsequent d-f^ are uniquely defined and can be numeri- 
cally generated with help of Eqs. ([^ through (^5)1 . At this time we shall put aside the discussion 
of the complex question of summability of the resulting multiple series in k^ and m? . 



5 The bubble-chain revisited 

Since the graph in Fig. 4 is of the order of 0(iV°), it appears in the higher-order diagrams only as 
one of the loops of the arbitrarily long 'bubble chain' shown in Fig. 5. Typically, the bubble chain 
is summed as a Dyson sum, even though (for qb, > 0) this generates a Landau-type singularity in 
the composed $2 field. The presence of this singularity in the j^{^)4 model was attributed to the 
expansion around a ground state that is unstable due to the effect of symmetry restoration |15) . 



12 




With massive fields (and also infinite-polynomial interactions) the pole still appears as a conse- 
quence of Borel nonsummability of the underlying series. The effects of Borel nonsummability were 
investigated in this context in and in jTS] based on the conjecture that that the Dyson sum 
is only a representative of a class of possible solutions which differ from each other by a function 
whose expansion in powers of the coupling constant is identically zero. 

Ref. [19 contains a prescription for finding a tachyon-free solution that differs from the (local 
theory's) Dyson sum by an addition of a nonlocal term which cancels the residue of the pole at fco, 
but in the perturbative expansion identically vanishes term by term. Since our present objective 
is restricted to the removal of divergences, it suffices to summarize that it is possible to construct 
a non-local effective model that is free of the Landau pole and which in the NLO, in the limit of 
p oo, behaves similarly to the jf{^)4 model. Furthermore, since (with the sole exception of the 
bubble-chain), there is no further infinite summations in higher orders in 1/A^, the introduction 
of a correction that eliminates the Landau pole could be regarded as an initial adjustment of the 
original Lagrangian in order to meet the requirement of consistency. 

An interesting property of the summed bubble-chain is that it contains a hidden subtraction 
which can help with the renormalization procedure. It can be best explained on the example of 
the simpler ■^{^"■^"y model. 

If we denote the contribution from a single bubble by I{k) (with all constants absorbed with 
the only exception of an overall factor of the coupling g) , and include in the sum the simple four- leg 
vertex, the result of the summation is 

1-g I(k) 

however, when the simple vertex is not included in the summation, the same result can be written 
as 



Since I{k) grows logarithmically with /c, in the high momentum regime the right-hand part of 
tends to —g. The single difference between the asymptotic behavior of graphs that contain the 
incomplete bubble chain and these that include only the simple four-point vertex is in the overall 
sign factor. This way the graphs containing the incomplete bubble-chain produce a subtraction 
for their simple- vertex counterparts. This kind of subtraction is not enough to reduce the number 
of divergences, but makes them milder by turning the overall logarithmic divergence of the graph 
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into a double-logarithmic one. For this reason the standard renormaUzation procedure still must 
be implemented. 

However, the hidden subtraction has an important beneficial consequence for the consistency of 
the dimensional renormalization which can handle only harmless poles whose residues are polyno- 
mials in masses and momenta, but not their logarithms [20j . In multi-loop diagrams the elimination 
of logarithmic poles is possible through mutual cancellations between graphs. In the early years 
of the dimensional renormalization, Collins [H] verified that in the model the cancellations of 
logarithmic residues indeed take place on the loop level. 

In contradistinction, in models there is no one-to-one correspondence between the number 
of loops and the order in powers of 1/A^. Therefore the loop-level cancellation argument of Ref. 
pi] no longer applies. However, the decomposition indicates that the logarithmic poles still 
disappear through the mechanism of implicit internal subtractions. Every effective four-leg vertex 
of every graph is made of a simple vertex and a bubble chain which are identical from the viewpoint 
of the renormalization procedure since their asymptotic forms differ only by an overall sign factor. 
Harmless poles can be removed by multiplicative renormalizations which can be evaluated from 
graphs including only the simple vertices. The logarithmic poles remain, but they cancel due to 
the mentioned sign difference. This effect is distinct, but consistent with the mechanism of mutual 
cancellations described in Ref. [3T]. 



6 Two-point function in NLO 

There are two types of NLO two-point graphs, shown in Fig. 6, that need consideration. As just 
explained, all we need to consider are graphs build with only simple four-leg vertices. When the 
same renormalization constants are applied to the graphs including the incomplete bubble chain, 
they cancel their singularities as well. 




Figure 6: NLO contribution to the two-point function. The second pair of graphs, shown in (b), 
is only an example; additional contributions are generated by inserting a bubble chain into one of 
the loops in Fig. [31 



Assuming that the asymptotic behavior of the propagator is canonical, in the zero order in 
combined powers of rri^ and in the vertices, the relevant (from the renormalization point of 
view) part of the two-loop self-energy diagram of Fig. 6a equals 
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The integral over the Feynman parameters in ([55]) contains a pole singularity produced by the 
divergence of the integrand at equal values of any pair of the parameters involved. The harmless 
poles of psp are removed by renormalization constants with 
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bi2 



1 (4^o¥ 

2N (47r)4 ' 
N (47r)4 ' 



(36) 
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The logarithmic poles cancel by the same mechanism as above. 

In the first order in combined powers of mass and momentum one encounters three types 
of vertex configurations, one with a simple multiplier of m|j.i|j, and two with factors including 
products of momenta. The extra actors of p originate either from an external leg and an internal 
propagator, or from a pair of internal propagators. The numerical value of the (harmless part of 
the) pole part of the diagram with a simple mass multiplier equals 



„2 m(2).3,(2) 

yR t2 l"oo ) "oi 
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The contribution from the diagram in which the momentum factor is produced by the external 
legs equals 



/^(2)x3j(2) 



1 



M -yi?-o ^ {2Tr)^^ J {p^-ml)[{q + k)^^ml][{p + qy -ml] 
It can be evaluated by an adaptation of the method of Ref. [21] which yields 



(38) 
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The integrand in pop is singular at a + /3 = and a + 7 = 0, but the transformation 



(39) 



a = px, (3 = p{l~x), 7=1- p. 



(40) 



converts the singularity at a + (3 — into a singularity at p = 0. In this limit the singular part of 
(|39|) equals 

(mg)"-3 



(m2)"-3 / dpdx(l -a:)p^-5 
^0 



n-A ' 



(41) 
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which is the pole part at a + /3 = 0. The pole part at a + 7 = is identical, and at /9 + 7 = the 
integral is regular; the pole exists, but at n = 6. 

The divergence is regularized by subtracting 



1 



1 



+ /3 a + 7 



(42) 



from the integrand. The integration can be performed exactly. The result is somewhat convoluted, 
but the numerical value of the divergent part of ((38)) equals 
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The last contribution comes from a graph in which the vertex multiplier involves momenta 
carried by two internal lines over which the field summation is carried out. The corresponding 
integral reads 
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The integration of the part that contains p'^Qp, in the numerator can be performed along the 
same lines. Its divergent part equals 
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The case of the part of the integral that contains p^ in the numerator is different since the 
integration produces not only -^^^ poles, but also an integral that behaves as dpp^^^ and is 
singular in the ordinary sense. However, the calculation can be easily carried out with use of the 
decomposition 



(p2 - rril)[{q + k)^ - ml] [{p + q)^ - ml] 
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The resulting pole part equals 
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Combining the divergent parts ([571) . I143p . (|1S1) . and (|¥71) . and using the LO information that 
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2d{Q' we obtain 
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The poles of need to cancel with analogous terms from the original Lagrangian (corrected 
with Zg, Zm, and Z$) and contributions from the LO diagrams with vertex coefficients dj./ so 
modified that they generate 0{1/N) corrections. So far any restrictions on relations between 
coefficients d'f^ were valid in 0((i)*') order only, and can be adjusted order by order in higher 



powers of -i-. In particular, in the NLO 



"■jk - "'jk + j^^jk ' 
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with the sole exception of c^qq which, as explained in Sec. 4, is exact, hence (5, 



(2) 
00 



0. In conse- 



quence, the LO renormalization constants remain unmodified. 

The residues of the five poles in (^5]) need to cancel with the residues of the poles generated 



by graphs represented by the first diagram of Fig. 2, with Sjj, included in in their vertices 



Additionally, conditions relating coefficients d'l' , i ^ 2 might be imposed as necessary. Extreme 



--jk 



(2) . 

care must be exercised when imposing additional relations between the d^j. coefficients, because 
LO constraints already have been imposed on them in Sec. 4. 

At this point we need to consider the NLO 0{g\) contributions from Z$ (i^Q^fc'^LQ, Z$ (f^^k'^mlL'^, 
and Z^ (i^Q^mpLQ. The case of the term is special since the LO two-point diagrams do not con- 
tribute to the fc^ terms. In consequence, the value of the coefficient djo^ in the nonlocal propagator 
must be such that the expression 
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is finite. Since (igo = 4, the pole cancels, provided that 
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While there is no restrictive conditions on ^ the coefficients 5\\ 



5(1) 
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c(2) 
'01 



and can be 



determined from the condition of cancellation of the pole parts and expressed as combinations of 
the LO coefficients of the Lagrangian. A straightforward calculation shows that 



5^^^ = 4.505d^V + 4.00l4o^ = 4.5054V " ll-9084o > 

= 1.6344V - 5.609dVV - 3.383dfo^ = 1.634^^ - 5.609dVV + 10.0684o\ 
4V = 0.8004V - 1-0004V + 16.3754o^ = 0.8004V - l.OOO^V " 48.7354o\ 

= -4.73l4V - 8.0094o = -4.73l4V + 23.8364V- (52) 
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In the limited context of this work there is no need for a detailed discussion of the contributions 
generated by the diagram in Fig. 6b. Its internal lines do not carry external momenta thus the 
integrations produce a series in powers of mass only. Mass renormalization can be approached 
in a manner similar to that used in the case of single- loop contribution from graphs in Fig. 2. 
Graphs with all vertex factors of d^^Q define the 0{g^) contribution to Zm- Further, in higher 
orders of p and m, the effect of the insertion of a three-loop 'blob' on a propagator line is the 

(2) 

same, independently of the values of k and / in the external dj^/ factor of the connecting vertex. 
In consequence, just as in the case of the one-loop LO contributions, it is possible to remove the 
divergences of all members of this class of self-energy graphs with a single multiplicative mass 
renormalization constant. 

Another amplitude that will be omitted for sake of brevity is the NLO six-point function 
which includes three types of graphs. One is essentially a bubble-chain with one of the four-point 
vertices replaced by a six-point vertex with extra legs being external. However, because of the 
identical topology, the cancellations of the related poles can be produced by a mechanism that 
closely parallels that which applies to the LO ($^)^ vertex. The other graphs are single-loop, with 
three ($^)^ or one ($^)^ and one (<I>^)'* vertex that produce six external legs. The former type 
of a diagram is finite in the local theory, but here it generates pole singularities due to additional 
powers of momenta in the vertices. The residua of all poles can be matched with the expansion 
coefficients of the ($^)'^ vertex so that the respective singularities cancel. 

7 Extension to higher orders, conclusions, and speculations 

Our analysis of the model based on the Lagrangian ([1]) demonstrates its viability as a candidate 
for an effective model of interaction of multiple scalars. The interaction has two components: 
a contact term (described by the local sector of £), and the nonlocal part of C, responsible for 
interactions mediated by the background parton environment. While the model was discussed in 
the context of interactions of colored scalars, other applications are also possible, e. g. scalar 
bound states in meson sector. Interestingly, the coupling of such hadronic molecules to mesons 
can also be modeled by a nonlocal interaction [52]. 

We did not attempt to discuss the renormalization properties of the amplitudes beyond the 
NLO of the expansion, and postponed the investigation of the behavior of the complete series 
of expansion in powers of momenta. While the number of free parameters in the model is sufficient 
to perform the presented procedures, and general patterns revealed in the LO and NLO in 
leads one to believe that they remain valid also in higher orders, one must be careful not to make 
overly optimistic claims. 

For example, the assumed property of the bare Lagrangian was that the asymptotic momentum 
behavior of the propagator is canonical. For this to be the case, it suffices that in the asymptotic 
region the nonlocal part of the propagator £5, given by Eq. ^ is suppressed by its canonical part, 
^Cs{p^ jTii?) — > (p^ — > 0), and :^Cs{p'^, m?) — > {w? — > 0). There is a danger of inconsistency 
though, since the requirement of ability to remove all types of divergencies with help of a single 
Z„i imposes relations between d^J.-* and d^^.'' . Some of the restrictions, like (PTjl can be regarded as 
only mild conditions of regularity, but others, like those discussed in Sec. 4 are more restrictive. 
While it might be possible to generate a set of jj^^ that would permit cancellation of the poles, 
maintaing the desired asymptotic behavior of the propagator Cs through higher orders in-i-is an 
entirely different question. 

Technically, one could maintain the canonical asymptotics by using unitary transformations to 
transform the violating terms to the vertices and repeating the procedure order by order inl/iV. 
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This however would not resolve any of the formal problems like global unitarity, which would 
require resummation of the entire series in and . There is a chance that such resummation 
might (at least in principle) be possible. The reason is that for the higher-order vertex to cancel the 
residues produced by the graphs that include a product of several lower-order vertices, its single 
coefficient must be equal but opposite to the sum of product of the lower order coefficients. The 
fact that the coefficients are of opposite sign suggests that the signs of the coefficients alternate 
(since the sums of products involved, this not necessarily has to happen order by order). Further, 
the magnitudes of higher coefficients compare with the products of lower order coefficients. If these 
are small, the expansion could possibly be summable within a finite radius of convergence. 

Another essential point is that consistency of the proposed scheme depends on the absence of 
logarithmic poles. The argument for their cancellation entirely relies on the presence of bubble 
chains in diagrams. Even though in higher orders in it is possible to construct graphs that do 
not contain a quartic vertex, and hence also no bubble chain exactly as discussed in Sec. 5, the 
procedure still works. The bubble chain is of the order 0{{j^)^) hence, following the prescription 
illustrated in Fig. 6, it can be attached to every <I)°<I)° propagator pair in higher-order vertices. 

At this juncture we are not entitled to make a claim that we have explained a mechanism of 
indirect interaction of scalars in a dense medium. Our goal was to investigate the structure of pole 
singularities. The virtue of the presented model is that it permits determination of finite constants 
in the Lagrangian from the assumption that pole parts of diagrams mutually cancel so that only 
the standard number of multiplicative renormalization constant is needed to yield finite results. 
The procedure presented here should be regarded only as a tool for generating an asymptotic series 
(in <1)°<1)°, p^, and m?) for an effective interaction, particularly applicable in the low-momentum 
transfer regime. 

Paradoxically, the inclusion of vertices with an unlimited number of external legs, that normally 
results in nonrenormalizability, was in fact essential for the renormalization procedure. Since the 
model was nonlocal, the presence of such vertices should not be interpreted as a dynamical feature 
of a multiparticle interaction, but rather as many-body correlations resulting from a number of 
intermediate interactions with the medium that eventually affect the scalars. In a nutshell, corre- 
lations between unstable bound states of particles of a renormalizable model should not render it 
nonrenormalizable. 

Acknowledgement: This work was supported in part by a grant from The Citadel Founda- 
tion. 



Appendix 

An alternative method of dealing with the divergences of the LO two-point amplitude is to 
proceed without renormalizing the nonlocality range, instead using a procedure applied to the 
four-point diagram. We opted for the first alternative for the reason of economy, and because it 
does not require any upfront constraints on the coefficients of the nonlocal part of the propagator. 
This makes it easy tomeet the requirement of the canonical asymptotic behavior of the LO propa- 

(2) 

gator. The alternative approach would require a prior definition of d^,; coefficients which are not 
determined until the next order of the 1 /N expansion. In the lowest order in , only I^i requires 
renormalization and the expression (jlOp correctly defines Zm- In the next order one renormalizes 
the sum of /qi -I- Iio . The part of /qi + Iio that is relevant for the renormalization procedure reads 
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^01 + /lo = r(i - ^) + (53) 
All graphs involved so far were of the first order in gj^, therefore the factor in (|53p can 

(2) 

be dropped. From the conventional point of view, separate renormalization constants for d-J are 
needed (which would imply nonrenormalizability in its traditional meaning), or imposition of the 
condition d^^-^ — — dxo' (which would imply triviality). However, the propagator ([2]) also contains 
the term 

lmldg^^'^{mlLl)^\ (54) 

and the cancellation of infinities is possible, provided there exists an appropriate relation between 
^02'' s-nd d^^-^ + d^o^. In our case this means that d^^^d^oo — d-m ^" ^lo^- Order by order, similar 
relations between d]^^ and d^^^^^_.^ also follow. However, in order to determine the exact form of 
the condition responsible for the cancellation of the pole of mass dimension (mg)^, we would need 
to consider the sum of all expressions ^ for which k + I = s. This would lead to a condition 
linking d^^^J to a weighted sum of c?^'^ which could not uniquely define d^^j^ . 
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